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1. Introduction 

The extragalactic background light ( EBL) is the background radiation field originating from 
starlight emission and its re-processing by interstellar medium integrated over cosmic history. The 
spectral energy distribution (SED) of the local EBL (at redshift z = 0) contains unique information 
about galaxy formation and evolution. Direct measurements of this SED are however very difficult 
because of foreground contamination. Strict lower limits have been derived from galaxy counts, 
and constraints on the EBL can also be derived using extragalactic very-high energy (VHE) y-ray 
sources. Eor a review see [2] and [3]. Indeed, VHE 7 rays interact with IR to UV photons via 
electron-positron pair production, resulting in an attenuated observed flux over fhe corresponding 
energy fhreshold. While fhis limifs fhe sfudy of exfragalacfic y-ray sources, if also provides a way 
fo probe fhe EBL over cosmic history. The aim is here to perform a measuremenf of fhe EBL 
densify based on y-ray specfra. The optical depfh for a VHE phofon of energy Ey fraveling fhrough 
a medium wifh EBL physical densify n{e,z) from a source af is: 

z{Ey,Zs)=cj^ ^yy (^y (1 ^ ^ M^ (1-1) 

where dt/dz= (^ 0 (1 + 2 )\/Dm (1 +z)^ , /X = 1 — cos ( 6 ) wifh 6 being fhe angle befween 

fhe phofon momenfa, ethr{Ey,z) = g ii(i+z) Befhe-Heifler cross secfion for pair pro- 

ducfion [1]. The absorption effecf on fhe infrinsic specfrum of an exfragalacfic source is expressed 
as ii>ohs{Ey) = EBL absorption fhen leaves a redshiff and energy dependenf im- 

prinf on fhe observed VHE specfra. To perform a measuremenf of fhe EBL SED, a significanf 
defection of fhis modulafion is necessary, which requires sfudying specfra wifh high signal fo noise 
rafio. 

Differenf models of fhe EBL have been developed (see e.g [4], [5], [ 6 ]) and fhey more or less 
agree on fhe overall local EBL shape : an opfical fo near infrared bump corresponding fo sfarlighf 
emission, a far infrared bump corresponding fo fhe fhermal re-emission of sfarlighf by infersfellar 
medium, in particular dusf, and a mid infrared "valley" befween fhese fwo bumps. The exacf prop¬ 
erties of fhe energy disfribufion in fhe mid infrared is however sfill nof well known and requires 
more observational inpuf. 

A previous EBL sfudy wifh H.E.S.S [7] was performed by simulfaneously tiffing fhe EBL 
opfical depfh and infrinsic specfra of exfragalacfic sources wifh a maximum likelihood mefhod. The 
shape of fhe EBL SED was fixed a priori choosing fhe model of [4], leaving only fhe normalizafion 
free. The overall fesf sfafisfics led fo a 8 . 8 -a defection of EBL absorption wifh respecf fo no 
absorption, wifh a normalizafion factor of 1.27 +q} 5 (sfaf) ±0.25(sys). 

The new analysis presenfed here follows a very differenf approach, focusing on fhe defermi- 
nafion of fhe shape of fhe EBL SED in addition fo ifs overall normalizafion. The sample of blazars 
used and fhe dafa analysis are described in Sec. 2. In Sec. 3, fhe mefhod fo access fhe shape of fhe 
EBL SED is presenfed. Einally, fhe resulfs of fhis preliminary sfudy are shown and discussed in 
Sec. 4. 
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2. H.E.S.S. data 

The high energy stereoscopic system (H.E.S.S.) is an array of five imaging atmospheric Cheren¬ 
kov telescopes located 1800 m above see level in the Khomas Highland, Namibia, detecting 7 rays 
ranging from ~ 100 GeV '. Due to the peaked pair-creation cross section around the energy thresh¬ 
old, this is precisely the corresponding sensitivity range for EBL absorption from ~ 0.3 /rm to 
< 100 /rm. In the present study, only data from the first phase of H.E.S.S. with four telescopes are 
considered. The sources used are bright H.E.S.S. blazars with a cut on the detection significance 
of 10 a : PKS 2005-489 {z = 0.071), PKS 2155-304 (z = 0.116), lES 0229+200 (z = 0.14), H 
2356-309 (z = 0.165), lES 1101-232 (z = 0.186), and lES 0347-121 (z = 0.188). This cut in sig¬ 
nificance ensures a large number of 7 rays from the source, a high signal to noise ratio and good 
reconstruction of energy spectra. The different data sets used are listed in Table 1. As in [7], for 
variable sources data are divided into smaller subsets restricted in flux range. 

The main analysis is performed using Model analysis [ 8 ], and the spectra were reconstructed 
using an unfolding technique [9]. The covariance matrix determined during the unfolding proce¬ 
dure is used in order to take into account the correlations between bins. A cross-check analysis 
is performed with ImPACT analysis [10] and an independent calibration chain, yielding consistent 
results. 


Data set 

z 

Emin-E,nax (TeV) 

lES 0347-121 

0.188 

0.20-6.5 

lES 1101-121 

0.186 

0 . 21 - 8.8 

H 2356-309 

0.165 

0.19-10.3 

lES 0229+200 

0.14 

0.69-10.2 

PKS 2155-304 (1) 

0.116 

0.23-10.1 

PKS 2155-304 (2) 

0.116 

0.21-9.6 

PKS 2155-304 (3) 

0.116 

0 . 21 - 12.0 

PKS 2155-304 (4) 

0.116 

0 . 21 - 11.1 

PKS 2155-304 (5) 

0.116 

0 . 21 - 6.0 

PKS 2155-304 ( 6 ) 

0.116 

0.26-15.0 

PKS 2155-304 (7) 

0.116 

0.26-13.9 

PKS 2155-304 ( 8 ) 

0.116 

0.21-12.9 

PKS 2005-489 (1) 

0.071 

0.26-8.9 

PKS 2005-489 (2) 

0.071 

0.26-11.1 


Table 1: Data sets used in this study with redshift and energy range covered by the observed 7 -ray spectra. 


3. Method 

3.1 Accessing the EBL shape 

The local, z = 0 EBL SED is described by splines constructed from a logarithmic grid in the 
(A, viv) plane in order to test a large variety of EBL shapes. Each spline is defined by knot points 

*For the fist phase of H.E.S.S.. This is brought down to a few tens of GeV with the fifth telescope. 
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on the grid as in [11, 12]. All tested shapes comply with the strict lower limits from galaxy counts 
(see references in [13, 3]), while the maximum is set relatively high over current upper limits in 
order to be conservative. Two shifted grids against each other are used to reduce the constraints on 
shapes imposed by the positions of the knots. Each shape is used to fit the data and is associated 
to a goodness-of-fit estimator in the form of The range covered by the grids covers the range 
of H.E.S.S. y-ray absorption sensitivity (~ 100 GeV - 20 TeV). In total 116,640 EBE shapes were 
tested. 



Figure 1: Grid setup to constmct EBL shapes. The black and grey dots are the knots for the two shifted 
grids. Also shown are lower limits from galaxy counts in green, upper limits from direct measurements in 
red extracted from [3], [13]. The dotted curve is the reference local EBL from the model of [4]. 

Eor each spectrum, a fit with EBE absorption is performed considering every shape i of the 
two grids, ^obs{Ey) = <!>,>,;The intrinsic (de-absorbed) VHE spectra are fitted with 
a power law and a log parabola, the best fit of which is selected considering the lower x^- Eurther 
improvements of the analysis will include consideration of other spectral shapes. The optical depth 
depends not only on the local z = 0 EBE density, but also on its evolution with redshift. This 
intrinsic evolution of the EBE is usually taken into account by introducing an additional scaling 
index to the cosmological photon density in the form of (1 with feyoi ranging from 1.2 

to 1.7 [12, 13]. Another approach is used here by taking the ratio of the density at redshift z to 
the density at redshift z = 0 from the model of [4]. It was checked that this does not introduce a 
strong model-dependence since a comparison with a fixed scaling index did nol yield significanl 
differences. Eurfher improvemenfs of fhe analysis will include differenl hypofheses. 

3.2 Combination of several spectra 

At first approximation, the EBE photon field should be essentially identical to all sources [14]. 
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Since each source provides a different combination of redshift and energy range, the results from 
the spectra of different sources can be combined to determine the EBL in a more significant way 
than using a single spectrum and over a wider range of energy, as a collective signature of EBL 
absorption. 

Eor each data set the envelope of all z = 0 EBL shapes inside the Xmin + ^ interval is retrieved. 
This envelope is considered as a 1 — a statistical contour of an actual measurement of the EBL 
SED. Note that this measurement is made at specific wavelengths defined by fhe grid spacing. 
The combinafion of fhe measuremenfs from each dafa sef is done as a combination of from 
independenf measuremenfs. This way, fhe combined cenfral value for each wavelengfh is fhe resulf 
of fhe average of fhe measuremenfs from each dafa sef weighfed by fhe corresponding 1 — a values. 

4. Results and discussion 

The resulf of fhis preliminary work fo access fhe EBL shape is summarized in Eig. 2. In fhaf 
figure, fhe grey area shows fhe combined besf-fif EBL SED shape from all considered dafa sefs, 
compared fo previous measuremenfs and consfrainfs, as well as fhe fiducial model of [4]. In order fo 
be conservafive fhe wavelengfh range presenfed here have been resfricfed fo mafch fhe pair creation 
fhreshold wifh fhe lowesf of fhe lasf poinfs in energy from our dafa sample, i.e 6.5 TeV. Considering 
fhe highesf poinfs of fhe dafa sample would lead fo a furfher coverage in fhe infrared. 
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Figure 2: Preliminary result of this work together with previous results, limits and the model of [4] 
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When comparing with previous measurements of the EBL density by H.E.S.S. it is important to 
note that the assumptions made are very different. Indeed in [7] the z = 0-shape and evolution of the 
fiducial model [4] was assumed and only the normalization was set free and fitted to the data. Here 
no prior is used on the considered z = 0 EBL shape and still the shape determined with our method 
is relatively close to the model of [4]. In Eig. 2, the dotted grey lines indicate the level of systematic 
uncertainties on the shape determination. Those have been obtained using different calibrations, 
analysis chains and spectral reconstructions. Part of the width of the systematic uncertainty comes 
from a jackknife analysis where some data sets were removed one by one. The final grey shaded 
area in Eig. 2 is drawn including all data sets. The main analysis leads to an EBL intensity around 
20 nW voT^ sr^ * (at the optical peak) slightly higher than previous measurements, at the level of 
upper limits. However, this excess is not significant given our current understanding of systematics. 

In the near future, improvements of the analysis will include consideration of other intrinsic 
spectral shapes. It is foreseen to include more sources at various redshifts and consider different 
hypotheses regarding the evolution of the EBL SED. Using spectra obtained with the second phase 
of H.E.S.S. with a lower energy threshold will also provide more leverage to constrain the short 
wavelength range of the EBL. 

5. Conclusions 

The results presented here are the first measurement of the local extragalactic background 
light SED shape with H.E.S.S., where previous studies dealt with its overall normalization. These 
preliminary results are based on first-phase H.E.S.S. data, and show that the sensitivity of the 
method allows to find results comparable to fiducial models in terms of shape, and broadens the 
previous H.E.S.S. measurements in terms of EBL wavelengths. 

This analysis is a first step towards a precise determination of the EBL SED with y-ray ab¬ 
sorption using H.E.S.S. data. This is a prerequisite for analyses aiming at the determination of the 
intrinsic VHE spectra of blazars and it is also closely related to alternative studies where the ef¬ 
fective mean free path of TeV 7 rays is modified (new fundamental physics like axions, secondary 
emission from cascades in the intergalactic medium, gravitational lensing,...). 
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